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Preface 


This  report  consists  of  five  chapters,  namely: 

I,1  Airfoils  in  Extreme  Ground  Effect,  pp.  2-6. 

2.  Automobile  Aerodynamics,  pp.  7-11. 

3.  Shallow-water  Waves  at  Discontinuities,  pp.  12-16. 

4.  Sliding  Sheets,  pp.  17-20. 

5.  Wings  Over  Water,y>p.  21-24. 


These  chapters  are  independent  of  each  other,  although  linked  by  the  “small-gap”  theme/  Chapters 
1-4  were  originally  presented  in  the  form  of  a  lecture  series  in  the  Department  of  Naval  Architecture 
and  Offshore  Engineering,  University  of  California  at  Berkeley,  during  February  and  March  1984. 
Chapter  5  formed  the  basis  for  a  presentation  to  the  panel  H5  (Analytical  Ship- Wave  Relations)  of 
the  Society  of  Naval  Architects  and  Marine  Engineers,  on  March  20th,  1984. 

Research  for  and  preparation  of  this  report  was  carried  out  during  a  period  spent  by  the 
author  in  the  Department  of  Mathematics,  Stanford  University,  January-May  1984,  sponsored  by 
the  Office  of  Naval  Research  under  a  contract  to  Professor  Joseph  B.  Keller.  Computing  support 
at  Stanford  was  provided  by  the  CLaSSiC  Project,  funded  by  the  Office  of  Naval  Research  under 
Contract  ONR  -82-K-0335  of  Professors  J.  Oliger  and  J.  Ferziger.  The  lecture  series  at  U.C. 
Berkeley  was  supported  by  the  Office  of  Naval  Research  Special  Program  for  Ship  Hydrodynamics, 
under  contract  N00014-84--K-0026,  to  Professor  Ronald  W.  Yeung. 

In  view  of  the  involvement  of  the  author  with  both  U.C.  Berkeley  and  Stanford  during  the 
time  of  its  preparation,  this  report  is  being  issued  in  identical  form,  under  separate  covers  of  both 


institutions. 


1.  Airfoils  in  Extreme  Ground  Effect 

In  this  chapter,  I  shall  concentrate  on  two-dimensional  irrotational  Sows  of  an  inviscid  incom¬ 
pressible  fluid  about  a  thin  airfoil,  as  in  Figure  1-1.  The  airfoil  is  very  close  to  an  impermeable 
plane  wall,  and  the  flow  is  in  the  first  instance  taken  to  be  steady. 

How  close  is  “eery  close ”  ?  Suppose  e  is  a  measure  both  of  the  airfoil’s  thinness  and  of  its  angle 
of  attack,  ie.,  if  L  is  its  length,  its  thickness  is  0(e).L,  and  its  angle  of  attack  is  0(e).  Then  we 
assume  that  the  clearance  or  gap  size  h  is  also  0(e).L.  That  is,  all  points  of  the  airfoil  surface,  top 
and  bottom,  are  not  more  than  0(e).L  from  the  wall. 


— —  v 


Figure  1-1.  Sketch  of  flow  situation. 

This  domain  of  closeness  is  indeed  “wry  c/ose”,  in  comparison  with  the  case  when  h  =  O(L), 
for  example.  The  aerodynamic  problem  of  h  =  0(£)  ground  effect  is  very  interesting,  and  deserves 
a  chapter  to  itself.  In  that  domain,  significant  (but  not  overwhelming)  flow  changes  occur,  relative 
to  the  case  of  an  airfoil  alone  in  an  infinite  fluid,  and  mathematical  and  computational  difficulties 
involving  numerical  solution  of  singular  integral  equations  have  to  be  overcome  before  solutions  can 
be  obtained  (see  eg.  Tuck  and  Newman,  1974).  We  say  no  more  about  that  case  here. 

In  the  case  of  present  interest,  h  =  0(e).L,  the  ground-effect  phenomenon  is  “  extreme  " ,  in  the 
sense  that  aerodynamic  effects  occur  due  to  the  presence  of  the  wall  that  are  much  greater  than 
those  that  occur  in  the  absence  of  the  wall;  specifically,  by  a  factor  0(e_1).  For  example,  the  lift 
on  a  thin  airfoil  at  0(e)  angle  of  attack  is  of  0(e),  whereas  the  lift  is  0(1)  in  the  extreme  ground 
effect  case. 

How  can  this  be,  when  every  facet  of  the  surface  of  the  airfoil  is  placed  at  only  a  small  0(e) 
angle  to  the  incident  stream?  If  the  wall  were  not  present,  such  a  body  would  cause  an  O(e) 
disturbance  to  the  stream,  leading  to  O(e)  pressure  changes  and  hence  the  usual  0(e)  lift  force. 

Indeed,  this  is  still  true  almost  everywhere,  and  leads  to  an  important  conclusion.  In  particular, 
the  flow  over  the  top  of  the  airfoil  is  still  an  0(e)  disturbance  to  a  uniform  stream,  and  does  lead 
to  O(e)  contributions  to  the  lift.  If  we  wished  to  solve  the  problem  to  an  0(e)  level  of  accuracy  (ie. 
to  find  0(1)  +  0(e)  contributions  to  lift,  consistently)  then  we  should  need  to  compute  the  0(e) 
disturbance  to  the  stream  over  the  top  of  the  airfoil. 

But  our  sights  are  set  somewhat  lower  (?).  We  are  here  content  with  the  0(1)  lift  contributions, 
and  hence  shall  neglect  all  0(e)  contributions.  This  has  a  remarkable  immediate  effect,  namely 
that  we  can  assume  that,  as  far  as  the  flow  above  the  airfoil  is  concerned,  the  airfoil  doesn’t  exist! 
That  is,  this  flow,  to  a  consistent  order  of  accuracy,  may  be  taken  as  an  undisturbed  uniform  stream 
of  magnitude  U.  Even  more  important  is  the  corollary  that  the  pressure  p  in  the  flow  above  the 
airfoil  is  unchanged  from  the  ambient  value  pa „  at  infinity. 
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So  much  for  the  flow  above  the  airfoil.  However,  some  fluid  may  pass  between  the  airfoil  and 
the  wall.  Unless  the  lower  surface  of  the  airfoil  is  exactly  parallel  to  the  wall  everywhere,  this  fluid 
will  be  moving  in  a  channel  of  variable  cross-section,  and  hence  its  velocity  will  be  non-uniform. 
Not  only  non-uniform,  in  fact,  but  in  general  varying  significantly  as  the  local  gap  width  changes. 
That  is,  we  must  expect  0(1)  flow  and  pressure  variations  in  the  gap,  and  this  is  what  leads  to  the 
0(1)  net  lift,  discussed  above. 

I  did  say  that  some  fluid  “may”  pass  beneath  the  airfoil,  not  “must”.  Indeed,  the  main  issue 
is  to  decide  just  how  much  fluid  does  so  pass.  The  answer  is  not  the  obvious  one  that  the  airfoil 
“scoops”  up  those  streamlines  that  would  have  passed  beneath  it  if  it  had  not  been  there.  What 
does  decide  this  issue? 

The  answer  is  that  trailing- edge  conditions  decide  it.  The  flow  must  detach  smoothly  from 
the  trailing  edge,  and  the  usual  Kutta  condition  of  aerodynamics  demands  that  the  pressure  be 
continuous  across  the  trailing  edge.  But,  given  our  consistent  0(1)  analysis,  with  the  conclusion 
that  the  above-airfoil  pressure  is  equal  to  ambient,  this  means  that  the  trailing-edge  value  of  the 
under-airfoil  pressure  must  also  be  ambient.  This  then  demands  that  the  fluid  velocity  in  the  gap 
must  return  to  the  free-stream  value  at  the  trailing  edge. 

Suppose  for  example,  that  the  airfoil  was  acting  as  a  simple  scoop,  as  above.  In  that  case, 
the  velocity  at  the  leading  edge  would  be  continuous,  and,  since  we  have  assumed  that  the  velocity 
takes  the  free-stream  value  outside  the  gap  region,  this  would  imply  that  the  leading-edge  value  of 
the  velocity  in  the  gap  was  also  of  the  free-stream  value.  This  fluid  would  then  change  its  velocity 
significantly  as  it  moved  along  the  gap,  and,  in  particular  would  emerge  at  the  trailing  edge  with 
a  velocity  that  is  in  general  not  equal  to  the  free-stream  value.  This  is  not  permitted;  hence  the 
scoop  concept  is  wrong,  and  we  must  allow  a  discontinuity  in  the  flow  velocity  at  the  leading 
edge. 

That  is,  somehow  the  leading  edge  knows  what  the  trailing  edge  is  doing,  and  selects  exactly 
the  right  amount  of  fluid  to  enter,  so  that  free-stream  conditions  will  prevail  at  the  trailing  edge. 
This  “selection”  process  involves  a  rapid  change  in  velocity  at  the  leading  edge,  since  the  velocity 
just  ahead  of  the  leading  edge  is  necessarily  (with  0(e)  error,  as  usual)  equal  to  the  free-stream 
value. 

This  situation  is  in  fact  not  unusual  in  aerodynamics,  just  more  extreme  in  its  consequences  in 
the  present  case.  The  Kutta  trailing-edge  requirement  always  performs  a  selection  process  among  a 
family  of  theoretically-feasible  flows.  In  unbounded  two-dimensional  aerodynamics,  selection  occurs 
via  a  choice  of  the  circulation  about  the  airfoil.  The  present  situation  is  no  different  in  that  respect; 
differing  strengths  of  the  leading-edge  velocity  discontinuity  represent  (indeed  are  proportional  to) 
differing  values  of  the  circulation  about  the  airfoil.  In  practice,  the  flow  must  adjust  itself  in  an 
unsteady  manner,  if,  for  example  it  is  started  from  rest,  until  the  Kutta  condition  is  satisfied  in 
steady  flow.  This  unsteady  development,  to  which  I  return  below,  involves  in  the  extreme  ground- 
effect  case,  a  change  from  an  initially  continuous  velocity  at  the  leading  edge  to  the  correct  value 
of  the  apparent  step  discontinuity. 

Of  course  there  is  no  actual  discontinuity.  If  we  were  to  examine  a  region  of  characteristic 
dimension  0{e).L  located  close  to  the  leading  edge,  we  should  find  a  rapid  but  not  instantaneous 
change  from  free  stream  to  a  different  uniform  velocity.  The  flow  in  this  leading-edge  neighbourhood 
is  indeed  very  interesting  and  readily  calculable  by  conformal  mapping  (see  eg.  Tuck  and  Bentwich 
1983),  since  the  boundaries  can  be  assumed  locally  plane.  This  local  flow  contains  the  necessary 
leading-edge  stagnation  point,  which  lies  at  a  distance  0(e). L  from  the  leading  edge,  on  the  airfoil’s 
surface  either  above  or  below  it,  as  in  Figure  1-2. 
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Figure  1-2.  Entrance- region  streamlines. 

Upstream  of  the  leading  edge  (locally)  there  is  an  apparent  source  which  models  rejection  of 
some  fluid  particles  that  otherwise  would  have  been  “scooped”  ;  the  source  is  negative  in  strength, 
ie.  becomes  a  sink,  if  an  excess  of  fluid  above  the  scoop  value  is  demanded.  Downstream  of  the 
leading  edge,  the  local  flow  asymptotes  to  a  uniform  stream,  of  exactly  the  right  magnitude  u(0) 
to  (eventually)  satisfy  the  Kutta  condition;  Figure  1-2  is  drawn  for  the  special  case  u(0)  =  U /2. 

The  action  then  takes  place  in  the  gap  itself.  What  action?  That  depends  on  the  geometry  of 
the  gap,  and  on  the  properties  of  the  fluid.  In  the  case  of  interest  here,  the  answer  is  very  simple, 
and  the  result  is  essentially  Venturi’s  law  for  one-dimensional  flow. 

That  is,  if  h  =  h(x)  is  the  clearance  or  gap  thickness  between  the  lower  surface  of  the  airfoil  and 
the  wall  at  station  x  measured  along  the  length  ,  then  the  flow  is  one-dimensional,  with  a  velocity 
u  —  u(x)  such  that  u (x).h(x)  is  constant.  Since  the  Kutta  condition  demands  that  u(L)  =  U,  this 
“constant”  is  determined  as  U.h(L)  ,  and  thus  u(x)  =  Uh(L)/h(x)  .  Hence  the  pressure  coefficient 
is  Cp(x)  ~  (p  -  Pao)/{?pU2)  —  1  -  h(L)7/h(x)3.  If,  as  discussed  above,  we  assume  p  =  p,*,.  i-e. 
Cp  =  0  ,  in  the  flow  above  the  airfoil,  the  net  lift  per  unit  span  is  just  F  =  \pU7  /  Cp{x)dx,  and 
the  centre  of  pressure  is  at  x  =  /  xCp(x)dxj  J  Cp(x)dx,  where  Cp(x)  is  given  above.  The  problem 
is  now  completely  solved,  given  h(x). 

It  is  of  interest  to  indicate  the  extremes  of  the  present  results.  Thus  Cp(x)  <  1  necessarily,  so 
F  <  \pU7L  ;  ie.  the  lift  coefficient  is  less  than  unity,  as  of  course  it  always  must  be  in  aerodynamics. 
However  (in  contrast  to  the  unbounded-fluid  case,  where  F  <<  L),  in  extreme  ground  effect, 
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a  lift  coefficient  close  to  unity  is  (at  least  in  theory)  achievable,  simply  by  letting  h(L)  vanish.  That 
is,  if  we  allow  the  trailing  edge  almost  to  touch  the  ground  (while  keeping  h(x)  still  significantly 
greater  than  zero  at  all  points  ahead  of  the  trailing  edge),  then  stagnation  conditions  will  apply 
over  most  of  the  gap,  and  hence  the  full  stagnation  pressure  will  be  felt  on  the  under-side  of  the 
airfoil.  Of  course,  the  velocity  must  return  to  the  free  stream  value  U  at  the  actual  trailing  edge  , 
but  the  net  flux  through  the  trailing-edge  gap  tends  to  zero  with  h(L). 

At  the  other  extreme,  an  arbitrarily-large  downward  force  (negative  lift)  is  theoretically  fea¬ 
sible,  by  allowing  h(x)  to  vanish  for  some  point  x  ahead  of  the  trailing  edge,  while  maintaining 
h(L)  >  0.  Again,  this  is  physically  obvious,  for  if  a  finite  flux  through  the  trailing  edge  is  required, 
and  if  h(z)  is  small  at  one  other  point,  then  a  large  velocity  is  needed  at  that  point,  which  gives  a 
large  negative  pressure  and  lift. 

These  considerations  and  others  are  illustrated  (see  Figure  1-3)  by  computing  the  net  lift  F 
and  centre  of  pressure  z  for  a  flat  plate  at  a  constant  angle  of  attack,  as  in  Tuck,  1981.  Various 
static  stability  results  can  be  deduced  from  the  behaviour  of  the  centre  of  pressure  as  a  function 
of  angle  of  attack.  Note  that,  for  small  angle  of  attack,  such  that  h(x)/h(L)  fa  1  everywhere  (ef. 
Widnall  and  Barrows  1970)  ,  one  can  linearise  the  theory  ,  and  it  is  then  easy  to  see  that  the  centre 
of  pressure  is  at  the  one-third  chord  point.  The  curves  of  Figure  1-3  for  a  —  1  are  relevant  to  the 
present  discussion;  those  for  a  values  other  than  unity  correspond  to  an  interesting  attempt  (Tuck 
1981)  to  account  for  trailing-edge  protuberances,  like  rudders.  The  parameter  0  used  in  Figure 
1-3  is  a  scaled  angle  of  attack;  the  related  parameter  A  is  the  contraction  ratio,  ie.  trailing 
clearance/leading  clearance. 


Figure  1-3.  Lift  and  center  of  pressure  for  flat  plates. 

Finally,  I  want  to  mention  some  extensions  of  the  ideas  discussed  here,  namely  to  unsteady  flows 
and  wakes,  and  to  non-rigid  walls.  Unsteadiness  is  of  importance  for  determining  development  of 
flows  from  rest,  for  (dynamic)  stability  considerations,  and  for  oscillatory  flows.  When  one  includes 
unsteadiness  in  the  extreme-ground-effect  theory,  one  must  also  concern  oneself  with  the  dynamics 
of  the  wake  vortex  sheet  shed  behind  the  airfoil  as  it  moves.  This  sheet  has  non-zero  strength  only 
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for  unsteady  flow.  The  question  (see  Tuck,  1978)  of  whether  such  a  non-trivial  vortex  sheet  will  or 
will  not  “roll  up” ,  is  an  interesting  one  that  can  be  answered  using  the  one-dimensional  theory. 

If  the  wall  is  not  plane,  but  is  still  fixed  and  impermeable,  the  theory  is  unchanged,  since  h(x) 
still  has  the  significance  of  the  net  clearance.  However,  an  interesting  generalisation  is  to  the  case 
of  a  moveable  wall,  in  which  there  may  be  a  second  coupling  between  the  pressure  p  in  the  air  flow 
in  the  gap  and  the  clearance  h.  This  type  of  extended  theory  is  discussed  further  in  Chapter  5. 
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2.  Automobile  Aerodynamics 

It  would  be  impossible  to  survey  the  held  of  automobile  aerodynamics  in  one  small  chapter  of 
this  report,  and  I  do  not  intend  to  attempt  this.  I  shall  concentrate  on  a  few  features  relevant  to 
the  main  theme,  and  hence  on  the  influence  of  the  ground  plane.  Recent  background  material  is 
contained  in  the  book  by  Sovran  et  at  (1978). 


Figure  2-1.  Flow  about  a  Ford  Capri. 

Figure  2-1  is  a  suitable  starting  point  for  discussion.  It  shows  the  results  of  a  two-dimensional 
irrotational  flow  computation  (Tuck,  1971)  on  an  automobile-like  profile.  In  fact,  this  profile  is  that 
of  the  Ford  Capri  of  that  year.  This  figure  has  been  reproduced  by  me  and  others  (e.g.  see  a  paper 
by  Bearman,  in  Sovran  et  01,1978)  subsequently,  and  seems  to  have  a  lot  to  say  about  the  general 
problem  of  automobile  aerodynamics.  That  is  a  somewhat  surprising  conclusion  at  first  thought, 
since  the  flow  about  a  real  automobile  is  neither  two-dimensional  nor  irrotational. 

Nevertheless,  the  computation  makes  those  assumptions,  and  uses  standard  (cf.  Hess  and 
Smith,  1964)  boundary-integral-equation  methods  to  solve  Laplace's  equation  exterior  to  an  arbi¬ 
trary  profile,  at  an  arbitrary  distance  from  a  plane  boundary.  In  fact,  two  distinct  boundary-value 
problems  are  solved.  In  the  first  there  is  assumed  to  be  a  uniform  stream  of  unit  magnitude 
at  infinity,  and  there  is  no  circulation  about  the  profile;  the  figures  written  outside  the  profile  are 
velocity  magnitudes  in  that  flow.  The  second  flow  (0j)  has  no  uniform  stream  at  infinity,  but  unit 
anti-clockwise  circulation  about  the  profile.  The  figures  inside  the  profile  are  velocity  magnitudes 
in  that  flow,  which  can  be  considered  to  have  been  generated  by  a  ficticious  unit  vortex  inside  the 
profile. 

The  general  solution  is  <j>  =  U<j>i  +/c^j,  where  U  is  the  vehicle  speed  and  k  the  actual  circulation 
about  the  profile.  Unfortunately,  we  don’t  know  the  value  of  k  in  advance,  for  an  arbitary  bluff 
body.  If  the  body  is  sufficiently  airfoil-like,  a  Kutta-like  condition  can  be  used  to  determine  k,  and 
we  return  to  that  case  later. 

In  the  meantime,  what  is  the  consequence  of  not  knowing  the  value  of  /c?  Interestingly,  there 
is  almost  no  consequence,  as  far  as  the  Sow  over  the  top  of  the  venicle  is  concerned,  as  can  be 
noted  from  the  relatively  small  inside  figures  over  the  top,  relative  to  the  outside  figures.  That 
is,  the  k  =  0  solution  is  adequate  to  describe  above-vehicle  flows.  It  is  tempting,  and  perhaps  even 
correct,  to  extend  this  conclusion,  and  to  claim  validity  for  a  zero-circulation  three-dimensional 
irrotational-flow  computation  (in  which  the  flow  beneath  the  vehicle  can  be  ignored),  as  a  predictor 
of  above- vehicle  velocities  and  pressures.  This  is  because  the  major  defect  of  such  a  computation  will 
be  lack  of  treatment  of  the  separated  wake,  but,  insofar  as  the  flow  ahead  of  the  wake  is  concerned, 
the  wake  serves  mainly  as  a  determinant  of  the  circulation,  and  we  have  seen  that  circulation  is  of 
minor  importance  over  the  top  of  the  vehicle. 
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In  contrast,  the  circulation  is  of  fundamental  importance  for  the  flow  beneath  the  profile,  as  is 
evidenced  from  the  relatively  large  inside  numbers  there.  Thus  wake  conditions  at  the  rear  end  are 
critical  in  determining  under-vehicle  flows,  and,  via  the  circulation,  the  net  lift  force  on  the  vehicle. 
In  three  dimensions,  induced  drag  will  also  be  profoundly  influenced  by  such  matters;  of  course,  in 
contrast  to  a  wing,  an  automobile  needs  no  lift,  and  hence  induced  drag  can  be  eliminated  if  the 
section-wise  circulation  is  zero.  The  pair  of  trailing  vortices  one  often  sees  behind  cars  on  a  dusty 
road  is  a  manifestation  of  this  induced  drag. 

This  point  is  illustrated  further  by  noting  that  since  4>  =  U<t>i  +  *4>2  is  a  linear  function  of  U 
and  k  ,  then  the  lift  F  per  unit  span  is  necessarily  a  quadratic  function  of  these  two  parameters, 
namely  F  =  FnU*  +  F^Uk  +  FttK*  ,  where  Fu ,  Fit,  Ftt  can  be  computed  by  suitably  integrating 
the  numbers  written  on  Figure  1.  For  the  Ford  Capri  profile,  it  turns  out  that  Fn  =  -O.lQpL7, 
F\ j  =  -1.27p,jPja  =  —0.3 Ap/L2.  Note  that  in  the  absence  of  ground  effect,  we  should  have 
necessarily  F\\  ~  Ft  a  =  0,  and  Fit  —  ~P  ,  by  the  Kutta-Joukowski  theorem.  The  term  in  Fn  is 
present  even  in  the  absence  of  circulation,  and  would  in  the  present  case  imply  a  small  downward 
force,  if  one  had  designed  the  vehicle  with  ic  =  0  to  eliminate  induced  drag. 

The  above  discussion,  while  restricted  to  two  dimensional  or  nearly  two-dimensional  flow,  has 
been  based  on  a  computation  that  does  not  assume  a  thin  body  or  a  small  gap.  The  qualitative  con¬ 
clusions  become  more  accurate  as  the  body  becomes  thin  and  the  gap  small,  however.  Nevertheless, 
so  long  as  the  body  is  bluff,  the  circulation  k  remains  unknown. 


r 


Figure  2.2.  Co-ordinate  sketch. 

Let  us  now  turn  to  the  case  of  thin  bodies  in  extreme  ground  effect  (Tuck,  1983).  This  is  the 
three-dimensional  generalisation  of  Chapter  1,  and  is  relevant  to  very  streamlined  and  low-to-the- 
ground  vehicles,  eg.  racing  cars.  Now  we  do  have  a  means  to  determine  the  circulation  k,  since 
a  smooth-detachment  or  Kutta  condition  must  apply  at  the  trailing  edge.  The  vehicle  is  assumed 
to  have  an  outline  curve  (as  seen  from  above)  of  the  form  T  =  +  Tr,  with  a  transition  point  T 
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between  leading  edge  IV  and  trailing  edge  IV,  as  in  Figure  2-2. 

The  mathematical  problem  for  the  flow  beneath  the  vehicle  is  to  find  a  velocity  potential 
<t>(x,y)  satisfying 

b(*2)  *«(*£)- 

inside  T,  subject  to 

<fi  =  Ux  (2.1) 

on  T*  ,and 

4>l  +  =  U 5  (2.3) 

on  Ft-  The  justification  for  these  equations  can  be  found  in  Tuck  (1983)  and  Tuck  (1980).  In  fact, 
equation  (2.1)  is  just  the  usual  shallow-water  equation  for  a  steady  flow  with  negligible  vertical 
velocity  in  a  region  of  non-constant  ‘‘depth’'  h(x,  y),  equation  (2.2)  matches  that  flow  to  the  incident 
stream  at  the  leading  edge,  and  equation  (2.3)  is  the  Kutta  condition  matching  to  the  ambient 
pressure  at  the  trailing  edge. 

A  difficulty  is  that  the  transition  point  T  is  not  known  in  advance.  That  is,  we’d  better  be 
clear  what  we  actually  mean  by  a  “leading  or  a  “ trailing"  edge.  It  seems  to  me  that  this  question 
has  received  insufficient  attention  from  aerodynamicists.  However,  that  is  excusable,  since  in  open- 
air  linearised  lifting-surface  theory,  T  is  necessarily  located  at  the  extremities  of  the  wing,  ie.  the 
leading  edge  always  changes  into  a  trailing  edge  at  the  wing  tips.  I  believe  that  this  is  true  in  all 
linearised  flows,  but  fails  in  the  present  case  because  the  problem  is  non-linear,  in  the  sense  that  0 
is  not  everywhere  close  to  Ux. 

Clearly,  a  leading  edge  is  one  where  fluid  passes  from  outside  the  contour  of  the  wing  to 
inside  it,  and  vice  vena  for  a  trailing  edge.  However,  this  doesn’t  quite  settle  it,  since  there  are 
two  insides,  namely  inside  above  and  inside  below.  Let  us  then  (at  least  in  the  purely  geometric 
sense)  be  quite  specific.  We  define  the  mean  horizontal-plane  velocity  vector  q  on  T  to  be  the 
average  of  the  vector  horizontal  fluid  velocities  above  and  below  the  wing  there.  Then  this  part  of 
T  belongs  to  IV  if  q.n  <  0  and  to  IV  if  qn  >  0,  where  n  is  the  unit  outward  normal  to  T  in  the 
horizontal  plane. 

In  classical  open-air  linearised  aerodynamics,  q  =  Ui  ,  since  the  perturbation  velocities  above 
and  below  the  lifting  surface  are  equal  and  opposite.  Hence  the  transition  point  T  between  leading 
and  trailing  edge  must  occur  where  the  contour  T  is  locally  parallel  to  the  free  stream,  ie.  at  an 
extremity.  The  same  applies  even  in  ground  effect,  if  we  linearise  the  problem  in  the  manner  of 
Widnall  and  Barrows  (1970). 

However,  in  the  fully-non-linear  problem,  this  is  not  necessarily  the  case.  Note  that  in  extreme 
ground  effect,  the  flow  beneath  the  vehicle  is  governed  by  (2.1)- (2.3),  but  the  flow  above  it  is  (to 
leading  order)  just  the  undisturbed  stream  <f>  =  l/x.  Hence,  in  carrying  out  the  averaging  process 
to  find  q  ,  we  use  17i  for  the  above- vehicle  component,  but  V<£  for  the  under- vehicle  component, 
and  since  <p  is  not  a  small  perturbation  of  U x,  we  no  longer  get  q  =  Ui.  In  fact  (cf.  Tuck  1983, 
Newman  1982),  the  local  streamline  at  T  necessarily  makes  an  angle  with  the  x-axis  that  is  twice 
the  angle  ie  a  by  T  . 

In  any  case,  it  is  not  difficult  to  build  determination  of  the  location  of  T  into  any  computa¬ 
tional  procedure  for  solving  the  boundary-value  problem  (2.1)-(2.3).  Since  the  problem  is  non-linear 
(because  of  (2.3)),  we  shall  be  forced  to  use  an  iterative  procedure,  and  we  can  simply  adjust  the 
location  of  T  at  each  iteration  step,  so  that  q.n  =  0  at  T.  Figures  2-3  and  2-4  show  sample 
computations  from  Tuck  (1983)  for  a  circular  lifting  surface,  at  positive  and  negative  angle  of 
attack,  respectively.  These  computations  were  done  with  a  special  (exponential)  variation  in  clear¬ 
ance,  such  that  (2.1)  possesses  a  (modified  Bessel  function)  family  of  eigensolutions,  which  can  be 


9 


used  to  simplify  the  task.  Work  is  proceeding  on  a  direct  attack  on  the  problem  for  more-realistic 
contours  and  clearances. 


Figure  2-3.  Circular  wing  at  positive  angle  of  attack. 

Meanwhile,  Figures  2-3  and  2-4  at  least  illustrate  some  important  phenomena  in  three- 
dimensional  extreme  ground  effect.  Figure  2-3  ( positive  angle  of  attack)  has  positive  pressure 
and  lift,  and  shows  fluid  being  deflected  outward  by  this  positive  pressure,  while  the  transition 
point  T  is  moved  forward  of  the  lateral  extremity  of  the  circle.  Conversely,  Figure  2-4  (negative 
angle  of  attack)  has  negative  pressure  (except  for  a  region  of  very  small  positive  pressure  near  the 
rear)  which  sucks  fluid  in  from  the  side,  and  the  transition  point  T  is  aft  of  the  lateral  extremity. 
Some  recent  racing  car  designs  have  exploited  the  ground  effect  for  negative  lift,  and  a  feature  of 
these  designs  is  the  use  of  “skirts*  or  lateral  barriers  near  the  extremities  to  prevent  the  inward 
flow. 


Figure  2-4.  Circular  wing  at  negative  angle  of  attack. 
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3.  Shallow*  Water  Waves  at  Discontinuities 


Shallow-water  waves  represent  perhaps  the  prototype  small-gap  flow,  the  water  depth  being 
the  “ small  gap ” .  The  water  is  shallow  if  h  is  small  compared  with  all  horizontal  scales  of  interest, 
of  which  perhaps  the  most  important  is  the  wavelength  A.  For  that  reason,  shallow-water  theory 
is  sometimes  referred  to  as  a long- wave”  theory,  a  terminology  that  I  have  long  opposed,  as  placing 
emphasis  on  too-narrow  an  interpretation  of  the  shallowness  property. 

In  the  present  lecture,  I  shall  only  discuss  linear  shallow-water  waves,  described  by  the  wave 
equation 


2b  -l.  3-(h2£)  -l  —lh21\  -  il 

dt  dx  dx  dy'dy  g  dt 3 


(3.1) 


where  the  depth  h  =  h(x,y,t)  in  general  can  depend  on  both  horizontal  space  co-ordinates  and 
time,  and  <f>  =  4>(x,  y,  t)  is  a  horizontal  velocity  potential.  One  can  interpret  cj)(z,  y,  t)  as  a  depth- 
averaged  three-dimensional  velocity  potential,  or  as  the  value  of  such  a  potential  at  the  bottom,  or 
at  the  free  surface,  etc.  Although  some  authors  make  a  big  deal  of  this  sort  of  distinction,  it  is  really 
rather  irrelevant,  the  point  being  that,  consistent  with  the  shallowness  assumption,  <j>  depends  only 
weakly  on  the  vertical  co-ordinate  z.  Note  that  the  actual  linearised  wave  elevation  is  —<j>t/g- 

The  bottom-motion  term  dh/dt  is  not  often  included  in  (3.1).  However,  an  important  appli¬ 
cation  where  it  is  of  crucial  importance  is  for  the  study  of  generation  of  tsunamis  by  earthquakes, 
with  h(x,y,  t)  as  the  given ( ?!)  ground  displacement  due  to  the  earthquake.  Of  course,  as  a  practi¬ 
cal  matter,  h(x,  y,  t )  is  known  only  sketchily  at  best.  In  any  case,  I  have  made  a  plea  (Tuck  1979) 
that  (3.1)  be  taken  seriously  as  the  governing  equation  for  tsunami  generation,  propagation,  and 
reception.  Effects  of  dispersion  and  non-linearity  are  of  very  little  significance  for  much  of  the 
history  of  these  extremely  long  and  low-amplitude  waves. 

Notice  that  if  there  is  no  time  dependence,  (3.1)  is  just  the  same  equation  as  used  in  Chapter  2 
to  describe  under- vehicle  air  flows.  The  famous  derivation  by  Friedrichs  (1948)  makes  it  clear  why 
a  similar  equation  will  always  hold  for  irrotational  incompressible  flows  in  a  small  gap.  Although, 
with  or  without  time  dependence,  (3.1)  is  not  entirely  straightforward  to  solve,  analytically  or  nu¬ 
merically  (see  again  Tuck  1979.  for  warnings  about  numerical  artifacts),  it  is  at  least  a  conventional 
wave  equation,  and  there  are  few  mysteries  once  we  accept  its  validity. 

However,  one  often  wishes  to  use  it  where  it  is  not  formally  valid,  ie.  where  there  are  horizontal 
length  scales  comparable  to  h.  In  particular,  it  is  not  entirely  obvious  what  are  the  correct  boundary 
conditions  in  some  cases,  especially  where  the  boundaries  are  not  vertical.  Stoker  (1957,p.420) 
gives  some  guidance,  but  never  addresses  the  difficulty  that  (3.1)  seems  to  lose  its  validity  near  the 
boundary. 

Perhaps  the  simplest  example  where  this  difficulty  arises  is  Lamb’s  (1932,  p.262)  treatment  of 
reflection  and  transmission  by  a  step  change  in  depth.  Suppose 


h  = 


*  <  0 
x  >  0 


and 


+  x  <  0 
Te*k+x-xot  x  >  Qi 


where 

gh±k$.  =  a2 

Physically,  this  corresponds  to  a  unit-amplitude  wave  in  x  <  0,  being  reflected  (co-efficient  p)  and 
transmitted  (co-efficient  r)  into  x  >  0.  Equation  (3.1)  is  satisfied  for  all  values  of  p  and  r,  and  our 
task  is  to  match  the  two  parts  of  the  solution  across  *  =  0  to  find  p  and  r. 


Lamb  uses  continuity  of  4>  and  of  h<f>Xl  namely 


1  +  p  =  t 


and 

h_k_(l  -  p)  =  h+k+r 


to  establish 


and 


P  = 


i-m* 

1  +  pi 


2 

1  + 


where  p  =  h+/h_  is  the  depth  ratio.  The  Lamb  continuity  conditions  are  similar  to  those  of 
Stoker.  Physically,  continuity  of  4>  implies  continuity  of  wave  elevation,  and  continuity  of  tuj>x 
implies  continuity  of  mass,  so  these  conditions  appear  rather  natural  and  correct. 

However,  continuity  of  h<j>x,  while  h  itself  jumps,  implies  that  jumps;  ie.  there  is  a  step 
discontinuity  in  the  horizontal  velocity,  which  seems  neither  natural  nor  correct  at  first  sight.  In 
fact  we  saw  this  difficulty  already  in  Chapter  1,  and  it  is  characteristic  of  situations  where  the 
small-gap  assumption  breaks  down  locally.  Here  the  breakdown  is  very  acute,  since  the  length 
scale  for  horizontal  changes  is  not  only  not  large  compared  to  h,  but  actually  vanishes  at  x  =  0. 

The  resolution  of  this  paradox  is  not  difficult  to  achieve  by  matching.  That  is,  we  abandon 
(3.1)  near  x  —  0,  and  seek  instead  a  solution  of  the  full  equations  of  motion,  which  means  in  the 
present  case  the  two-dimensional  Laplace  equation  (there  being  no  y-dependence)  in  the  vertical 
plane,  namely 

^s*  +  =  0 


for  i  =  0(h),z  =  O(h).  But,  in  such  a  spatially-stretched  region,  the  free  surface  appears  locally 
rigid  (ie.  <t>,  =  0  on  z  =  0)  a  conclusion  that  can  be  derived  formally  from  the  exact  free-surface 
condition,  for  waves  with  A  >>  h,  ie.  wavelengths  far  in  excess  of  the  length  scales  of  interest 
locally.  The  local  problem  is  thus  a  simple  streaming  flow  through  a  rigid- walled  channel,  of  step¬ 
changing  cross-section,  as  in  Figure  3-1,  and  can  be  solved  easily,  eg.  by  conformal  mapping.  In 
fact,  this  qualitative  picture  applies  for  a  general  depth  change,  not  necessarily  restricted  to  a 
sudden  jump,  but  in  the  general  case  the  numerical  solution  may  not  be  so  easy. 


Figure  3-1.  Sketch  of  local  flow  near  depth  change. 
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The  result  (see  Tuck  1376)  is  that  we  confirm  Lamb’s  jump  conditions,  and  in  addition  obtain 
a  first  correction  to  the  ^-continuity  condition,  namely  that  there  is  a  (small,  0(h))  jump  in  the 
value  of  <p  across  x  =  0,  of  magnitude 

A  <f>  =  2  C.h4>x 

(recalling  that  h</>x  is  continuous),  where  C  is  a  uniquely-determined  co-efficient  depending  on  the 
detailed  shape  of  the  depth  change.  For  the  specific  case  of  a  step  as  in  Figure  3-1,  we  have 


2j rC  = 


log 


M-l 


-  2  log 


4(i 


-  i' 


The  effect  on  the  Lamb  problem  is  that  the  reflection  co-efficient  is  corrected  to 

p  =  i^iil  -  4 4CJ^-k+h+  +  0(h3)] 


Note  that,  (for  p^l),  the  new  0(h)  term  affects  only  the  phase  of  the  reflection  co-efficient;  the 
first  correction  to  its  amplitude  arrives  only  at  the  next  0(/»2)  correction.  Similar  but  less  general 
results  were  obtained  by  Bartholomeuss  (1958). 


Figure  3-2.  Sketch  of  flow  for  paddle  wavemaker. 

The  ideas  behind  the  above  matching  can  be  generalised  in  a  number  of  other  shallow-water 
problems.  For  example,  in  my  own  work,  I  have  studied  paddle-type  wave-makers  (Tuck  1974) 
that  do  not  occupy  the  whole  water  depth,  thereby  leaving  a  bottom  gap  as  they  oscillate.  As  in 
the  Lamb  problem,  one  can  match  a  solution  of  (3.1)  (in  this  case,  a  purely  out-going  wave)  to  a 
local  rigid-firee-surface  flow  near  the  paddle,  as  in  Figure  3.2.  This  leads  to  predictions  of  the  loss 
in  efficiency  of  the  paddle,  as  a  function  of  the  gap  size.  If  the  gap  size  varies  with  time  during 
the  stroke,  as  for  example  it  would  if  the  paddle  were  swinging  like  a  pendulum,  the  corollary  is  a 
distorted  wave  form,  in  spite  of  the  basic  linearity  of  the  problem. 

As  part  of  a  general  study  of  water-wave  (not  necessarily  shallow)  transmission  through  small 
holes  in  vertical  barriers,  the  shallow-water  special  case  was  treated  in  an  Appendix  to  Guiney  et 
al  (1972)  by  similar  methods.  The  transmission  co-efficient  that  results  can  be  written 

r  =  1/(1  +  ikC) 

where  k  =  2ir/A  is  the  wave  number,  and 

_  D  .  a  .  irh. 

C  =  -V!ogl’rD‘mD, 
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in  the  local  geometry  of  Figure  3-3. 


Figure  3-3.  Local  flow  through  a  hole  in  a  wall. 

A  recent  study  (Tuck  1980)  of  waves  incident  upon  a  basin  partially  protected  by  a  reef  as  in 
Figure  3-4,  used  similar  concepts.  The  results  for  the  natural  frequencies  of  the  basin  as  a  function 
of  gap  size  compared  well  with  model  tests  (Tuck  et  aL,  1980). 


y 


Figure  3-4.  Flow  geometry  for  submerged  reef. 

Finally,  I  should  point  out  that  a  slender  ship  moving  in  shallow  water  provides  a  good  example 
of  an  apparent  discontinuity  in  application  of  (3.1)  (or  its  equivalent  for  some  steady  flows).  Again, 
one  has  to  match  with  a  local  flow  near  each  cross-section  of  the  ship,  in  which  the  free  surface  is 
replaced  by  a  rigid  wall.  Such  problems  are  discussed  in  Tuck  1976  and  1978;  see  also  Yeung  and 
Tan  1980. 
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4.  Sliding  Sheets 


G.  I.  Taylor’s  famous  movie  (1967)  on  low  Reynolds’  number  fluid  dynamics  contains  an 
entertaining  demonstration  of  paper  sheets  sliding  easily  over  a  smooth  table  top,  an  experiment 
repeatable  by  anyone.  Taylor  described  this  phenomenon  as  an  example  of  fluid  lubrication,  which 
it  certainly  is.  What  is  far  from  obvious  is  that  the  viscosity  of  the  fluid  plays  any  role.  That  is, 
could  the  paper  be  ‘‘flying’’ ,  like  an  airfoil  in  ground  effect? 

I  used  this  example  often  when  I  first  worked  in  this  area,  and  was  originally  quite  confident 
that  the  answer  to  the  above  question  is  “Yes”.  If  so,  the  inviscid  analysis  (e.g.  from  Chapter  1, 
for  a  sheet  of  large  aspect  or  width/length  ratio)  predicts  that  if  h(x),0  <  x  <  L,  is  the  clearance, 
then  the  pressure  under  the  sheet  is  given  by 

P  =  Pao  +  \pUi[l-h(L)i/h(x)i ]  (4.1) 

from  which  we  can  obtain  the  lift  F  and  moment  M  per  unit  width.  However,  there  are  two 
difficulties: 

(i)  F  <  \pU‘1L  ,i.e.  there  is  an  upper  limit  on  the  weight  of  sheets  that  can  fly. 

(ii)  M  0  about  the  centre,  i.e.,  a  uniform  sheet  cannot  be  in  pitching  equilibrium. 


Figure  4-1.  Flexible  sheets,  highest-flying  bending  modes. 

The  latter  is  a  major  difficulty,  and  I-  attempted  to  resolve  it  within  the  inviscid  theory  by 
allowing  the  sheet  to  bend  (Tuck,  1982).  The  resulting  analysis  is  an  interesting  exercise  in  non¬ 
linear  coupled  aerodynamics  and  elasticity,  and  an  infinite  number  of  equilibrium  configurations 
was  obtained.  Figure  4-1  shows  the  first  (highest-flying)  mode,  for  various  weights.  The  scale  of  the 
figure  is  inversely  proportional  to  the  stiffness  D,  and  hence  this  theory  predicts  that  a  nearly-rigid 
sheet  must  fly  very  dose  to  the  ground. 

However,  air  is  viscous,  and  in  particular,  the  Reynolds’  number  based  on  the  gap  clearance  h 
is  not  necessarily  large.  If  viscosity  ft  is  dominant,  we  might  expect  a  fully-developed  Couette-type 
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flow,  with  a  nearly-linear  velocity  profile  across  the  gap,  modulated  by  corrections  that  vary  slowly 
along  the  gap.  This  is  typical  of  “lubrication  theory” ,  and  leads  (for  a  flat  plate  at  angle  of  attack 
o)  to  a  pressure  distribution  given  by 


P  =  Poo  + 


3 fiax(L  -  x) 


(4.2) 


This  result  is  only  true  for  sufficiently  small  Reynolds’  number,  in  contrast  to  (4.1),  which  is  only 
true  for  sufficiently  large  Reynolds’  numbers.  Equation  (4.2)  also  has  a  defect  like  (ii)  above,  in 
that  the  moment  about  the  center  point  is  non-zero  —  indeed  the  sign  is  opposite  to  that  predicted 
by  equation  (4.1). 

So  what  actually  happens  to  a  non- flexible  uniform  sheet?  It  cannot  be  in  pitching  equilibrium 
at  either  high  or  low  Reynolds’  number.  Fortunately,  the  sign  variation  in  moment  at  these  extremes 
suggests  that  there  will  exist  an  intermediate  Reynolds’  number  where  the  moment  about  the  center 
is  zero.  I  now  give  a  brief  summary  of  a  study  of  this  problem  by  Tuck  and  Bentwich  (1983). 

The  Reynolds’  number  of  interest  is  actually  a  special  one,  namely 


_  Uh(° )a 

vL 


(4.3) 


which  is  a  small  multiple  h(0)/L  of  the  Reynolds’  number  based  upon  clearance  h,  and  an  even 
smaller  multiple  ( h(0)/L )J  of  that  based  on  sheet  length  L.  Hence  if  R  =  0(1),  both  the  latter 
Reynolds’  numbers  are  large,  and  an  in  viscid  theory  applies  in  all  regions  of  flow  except  the  main 
gap.  In  particular,  the  entrance-region  flow  is  in  viscid,  and  in  fact  is  the  same  as  Figure  1-2. 

But  we  are  mainly  interested  in  what  is  happening  in  the  main  gap  region,  and  it  is  not  hard 
to  show  by  asymptotic  stretching  that  the  appropriate  equation  is  the  boundary-layer  equation 

~  V'xV'w  =  “  ~P'(r)  +  ^vvv  (4-4) 


for  the  stream  function  ij)(x,  y).  We  have  to  solve  equation  (4.4)  subject  to  no-slip  boundary 
conditions  on  both  walls;  this  gives  four  conditions  for  a  3rd-order  equation,  enabling  simultaneous 
determination  of  the  (unknown)  pressure  p(x). 

Before  we  can  start  solving  (4.4)  we  need  “initial”  conditions  at  the  starting  station  x  =  0,  and 
“exit”  conditions  at  the  station  x  =  L  where  the  gap  flow  merges  with  the  inviscid  outer  flow  above 
the  sheet.  The  former  is  a  requirement  that  the  flow  match  an  entrance  flow  as  in  Figure  1-2  ,  with 
an  arbitrary  velocity  u(0)  at  x  =  oo,0  <  y  <  h( 0),  i.e.,  the  entrance  profile  is  a  uniform  stream 
of  to-be-determined  magnitude.  How  is  u(0)  determined?  The  pressure  p(x)  must  commence  with 
the  Bernoulli  value 

P(0)  =  Poo  +  ^ Pu 2  ~  (4-5) 

at  x  =  0,  and  end  with  the  free- stream  value 


P(L)  =  Poo  (4-6) 

at  x  =  L.  Equation  (4.6)  is  the  “exit”  condition  that  finally  fixes  u(0). 

We  now  have  enough  conditions  to  solve  (4.4).  The  resulting  boundary-value  problem  is  an 
“inverse”  problem,  in  the  terminology  of  Keller  (1978),  since  p(x)  is  not  prescribed  as  it  is  in 
the  more  usual  external  or  “direct”  boundary-layer  problems.  However,  there  are  many  practical 
numerical  methods  available,  and  we  can  consider  this  problem  as  easily  solveable. 
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Figure  4-2  shows  a  typical  set  of  streamlines  for  a  case  with  h(0)/h(L)  =:  4  with  a  very  tight 
constriction.  In  fact  whenever  this  “contraction  ratio”  exceeds  2,  there  is  a  separated  “bubble” 
attached  to  the  upper  (sloping)  face,  within  which  there  is  reversed  flow.  Nevertheless,  the  compu- 


Figure  4-3  displays  the  main  pressure  results  for  the  flow  of  Figure  4-2.  These  figures  show  the 
inviscid  (4.1)  and  lubrication  (4.2)  approximations,  together  with  a  useful  “linearised”  approxima¬ 
tion  described  in  detail  in  Tuck  and  Bentwich  (1983).  The  curves  labelled  “exact”  are  numerical 
boundary-layer  computations  to  high  precision. 


o  o.:  o.4  o.6  o.8  i.o 

x/L 

Figure  4-3.  Pressure  results,  flow  of  Fig.4-2. 


Figure  4-4  is  a  final  result,  obtained  by  integrating  this  pressure,  for  those  cases  where  the 
centre  of  pressure  is  at  mid-chord.  This  figure  gives  the  “design  charts”  for  steady  flight  of  uniform 
sheets  in  pitching  equilibrium. 


Figure  4-4.  Lift  and  angle  of  attack  for  uniform  sheets. 
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5.  Wing  a  Over  Water 


Suppose  we  have  an  extreme  ground-effect  air  flow  about  a  thin  airfoil  like  that  of  Figure  5-1. 
In  contrast  to  the  flows  treated  in  previous  chapters,  the  ground  plane  is  neither  a  ground  nor  a 
plane.  Instead,  it  is  supposed  to  be  a  free  water  surface,  whose  main  property  is  deformability 
in  response  to  aerodynamic  pressure.  We  assume  now  that  h(x)  denotes  the  height  of  the  under¬ 
surface  of  the  airfoil  relative  to  the  undisturbed  water  surface,  and  rj(x)  is  the  elevation  of  the 
disturbed  water  surface. 


Figure  5-1.  Sketch  of  airfoil  near  water. 


A  one-dimensional  air  flow  is  assumed  to  occur  in  the  small-gap  region,  and  hence  the  pressure 
distribution  varies  like  the  negative  inverse  square  of  the  local  clearance,  namely  as 


P(x)  =  Poo  +  ^PaU2\  1  -  ( 


h(x)  -  rj(x)  ] 


(5.1) 


where  is  the  ambient  or  atmospheric  pressure  at  infinity,  and  pa  is  the  density  of  the  (incom¬ 
pressible)  air. 

We  need  one  more  equation  relating  p(x)  and  rj(x).  First  let  us  assume  (Tuck  1984)  that  this 
equation  is 

p(z)  =  Poo  -  Pwm(x),  (5.2) 

i.e.  the  pressure  is  hydrostatic,  with  g  as  the  acceleration  of  gravity  and  pw  as  water  density.  In 
effect,  we  are  assuming  that  the  water  does  not  move  at  all,  but  its  surface  merely  deforms  statically 
in  response  to  the  aerodynamic  pressure.  Some  direct  justification  of  this  type  of  assumption  is 
provided  in  Tuck  (1975). 

Before  looking  at  the  consequence  of  (5.2),  it  is  instructive  to  note  that  the  relationship  (5.2) 
is  a  simple  linear  representation  of  potential  energy  storage  by  the  deformable  surface  y  =  tj(x), 
and,  as  such,  could  be  achieved  in  other  ways.  For  example,  suppose  the  “water”  were  replaced 
by  a  type  of  elastic  medium,  known  in  the  theory  of  foundations  as  a  “Winkler  soil”  (Selvadurai 
1970),  consisting  of  a  distribution  of  independent  linear  springs.  Then,  again,  there  would  be  a 
linear  pressure-versus-displacement  relationship  like  (5.2),  with  a  “modulus  of  sub-grade  reaction” 
or  spring  constant  per  unit  area,  of  magnitude  pw9-  A  qualitatively  similar  result  also  holds  if  we 
allow  the  air  to  be  compressible,  even  if  the  ground  is  still  rigid  and  plane,  since  again  we  should 
be  providing  potential  energy  storage. 
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In  any  case,  if  we  accept  (5.2)  and  combine  it  with  (5.1),  we  can  eliminate  r?(i),  obtaining  a 
single  non-linear  algebraic  equation  for  the  pressure  p(x),  namely 

h{x)-[h(L)-Tj(L)}(l-P(x))-i  (5.3) 

*  PW  9 


where 


P(»)  ~  P 
(paU*  ' 


(5.4) 


Equation  (5.3)  can  be  re-written  as  a  cubic  polynomial  equation  for  P(x),  and  solved  easily 
for  any  given  h(x).  There  are  special  significances  to  be  attached  to  the  parameter  choices  at  which 
the  discriminant  of  this  cubic  vanishes,  i.e.  where  the  number  of  real  roots  of  (5.3)  changes  from  one 
to  three,  and  this  is  discussed  further  in  Tuck,  1984.  In  particular,  physically-acceptable  solutions 
exist  if  h(x)  >  h(L)  whenever  x  <  L  ,  i.e.  the  trailing  edge  is  lower  than  any  other  point,  while  at 
the  same  time  F  <  1,  where 


f> ,  SLe*. 

gh(L)  pw 


(5.5) 


The  parameter  F  is  a  kind  of  modified  clearance-based  Froude  number,  and  the  above  requirement 
states  that  “sub-critical”  flows  exist  at  positive  angle  of  attack.  Figure  5-2  shows  the  computed 
free  surface  shape  *7(2)  for  a  flat  plate  in  some  such  cases. 


x/L 


Figure  5-2.  Hydrostatic  free  surface  under  flat  plats. 

In  more  recent  work  (Grundy  and  Tuck,1984),  the  hydrostatic  law  (5.2)  has  been  replaced  by  a 
true  hydrodynamic  model  of  the  water  surface.  The  simplest  way  to  achieve  this  is  to  note  that  the 
right-hand  side  of  (5.3)  is  just  the  water  surface  elevation  q(x);  if  (5.2)  does  not  hold,  but  instead 
we  allow  hydrodynamic  effects,  we  simply  replace  the  current  right-hand  side  of  (5.3)  by  the  known 
q(x)  produced  by  a  given  pressure  distribution  p(x). 

Since  the  hydrodynamic  disturbance  is  small,  because  the  free- surf ace  slope  rj'(x)  is  everywhere 
small,  the  free- surface  boundary  condition  can  be  linearised  (for  the  flow  in  the  water,  if  not  in 
the  air),  and  the  solution  of  the  hydrodynamic  problem  is  as  given  for  example  in  Wehausen  and 
Laitone  (1960;  see  also  Tuck  1982).  Thus  (5.3)  becomes 

M*l  -  MO  - 1(01(1  -  ewr*  =  -  f rn,  (5.8) 
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where 


K'{ x)  =  -/'(*)/*  -  2smx£f(x)  (5.7) 

with  f'{x)  as  an  (even)  auxiliary  function  for  the  sine  and  cosine  integrals,  and  H  the  Heaviside 
step  function.  The  problem  is  now  reduced  to  a  non-linear  integral  equation  for  P(x). 

An  interesting  paradox  presents  itself.  We  might  hope  (even  expect?)  that  (5.6)  reduces  to 
(5.3)  as  U  — ►  0,  i.e.  as  hydrodynamic  effect  become  weaker.  This  almost  happens!  If  we  let  the 
length-based  Froude  number  Fi  =  U(gL)~$  tend  to  zero  in  (5.6),  the  right-hand  side  tends  to  that 
of  (5.3),  but  only  if  both  P(0)  and  P(L)  also  tend  to  zero  in  that  limit.  That  is,  the  pressure 
distribution  must  return  to  atmospheric  at  both  ends. 

In  fact,  in  solving  (5.6)  for  general  non-zero  Fi,  we  can  and  do  demand  already  that  P(L)  =  0, 
i.e.  we  enforce  a  Kutta  condition  that  the  pressure  return  to  atmospheric  at  the  trailing  end.  This 
is  a  freedom  available  to  us  since  the  end  clearance  h(L)  -  rj(L)  that  appears  explicitly  in  (5.6)  is 
an  extra  unknown,  and  in  practice  we  allow  its  value  to  vary  until  P(L)  =  0.  However,  we  have  no 
reason  to  expect  that  P(0)  =  0  at  any  Pl,  and  in  particular  there  appears  no  numerical  evidence 
that  (in  general)  P(0)  — *  0  as  Pj,  — ♦  0.  Hence  it  must  be  exceptional  for  solutions  of  (5.6)  to 
approach  solutions  of  (5.3)  as  Pi  -*  0. 

In  any  case,  efforts  are  underway  to  solve  (5.6)  numerically.  I  show  here  (Figure  5-3)  perhaps 
the  most  interesting  of  the  preliminary  outputs.  These  results  are  for  a  flat  plate  at  various  positive 
angles  of  attack  in  sub-critical  flow  at  P  «  0.7,  as  measured  by  (5.5).  Hence  the  hydrostatic  theory 
gives  results  like  those  of  Figure  5-2,  which  have  the  property  that  P(0)  ^  0.  Figure  5-3  now  gives 
computed  values  of  P(0)  as  a  function  of  Pl,  at  fixed  P.  There  is  no  tendency  for  P(0)  to  vanish 
as  Ft  — *  0,  nor  in  fact  for  P(x)  as  computed  by  this  program  to  approach  that  computed  from  the 
hydrostatic  theory,  at  any  x. 


Figurs  5-3.  Leading-edge  pressure  oa  s  Sat  plate. 


Another  remarkable  feature  of  the  finite-  Fi  results,  displayed  in  Figure  5-3,  is  that  they  are 
non-unique  in  some  parameter  ranges!  For  example,  at  about  P&  =  0.25,  a  vertical  line  intersects 
the  curve  labelled  h(0)/h(L)  =  1.2  in  three  places.  I  must  confess  that  we  are  presently  only  able 
to  compute  two  of  these  three  solutions,  the  dashed  portion  of  the  curve  being  inferred.  There  are 
obviously  a  number  of  things  happening  here  that  are  not  easy  to  understand! 
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On  a  more  positive  closing  note,  Figure  5-4  shows  some  results  that  do  approach  the  hydrostatic 
theory  as  Fl  -*  0.  These  are  for  a  (parabolic)  curved  airfoil,  which  is  fore-and-aft  symmetric.  In 
such  a  case,  the  hydrostatic  theory  demands  a  fore-and-aft  symmetric  pressure  distribution.  Thus, 
if  P{L)  =  0,  necessarily  P(0)  =  0  also,  at  Fl  =  0.  The  numerical  work  for  finite  FL  confirms  this; 
although  there  is  no  fore-and-aft  symmetry  of  P{x)  so  long  as  Fl  ^  0,  such  symmetry  is  approached 
as  Fl  — »  0,  and,  better  yet,  the  hydrostatic  results  are  approached  in  that  limit  for  all  z.  Note 
of  course  that  this  limit  is  still  a  somewhat  singular  one,  since  the  Sow  at  any  Fl  ^  0  contains 
waves,  and  these  waves  become  shorter  and  shorter  as  Fl  — *  0,  even  though  (in  the  present  case  of 
a  fore-and-aft  symmetric  curved  body,  if  not  for  a  fiat  plate)  their  amplitude  is  tending  to  zero. 


Figure  5-4.  Free  surface  below  symmetric  curved  body. 
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